Mimicking the probability distribution of a two-dimensional Grover walk with a 

single-qubit coin 
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The non-localized case of the spatial density probability of the two-dimensional Grover walk can be 
obtained using only a two-dimensional coin space and a quantum walk in alternate directions. This 
significantly reduces the resources necessary for its feasible experimental realization. We present a 
formal proof of this correspondence and analyze the behavior of the coin-position entanglement as 
well as the x-y spatial entanglement in our scheme with respect to the Grover one. Our scheme 
allows us to entangle the two orthogonal directions of the walk more efficiently. 



The concept of a random walk, as the mathematical 
description of a sequence of random steps, has a large 
number of applications in different fields, ranging from 
physics to computer science and from economics to biol- 
ogy Its analogue in quantum theory is the so-called 
quantum walk which, due to the unitarity requirement, 
brings with it certain changes 0, H|: while in a clas- 
sical random walk the position of the particle (walker) 
is described by a probability distribution, in a quantum 
walk the particle is in a superposition of quantum states. 
It has been shown that any quantum circuit can be effi- 
ciently simulated using quantum walks Q and interesting 
proposals to implement quantum walks in different physi- 
cal systems have been put forward 0, @ . An equivalence 
between the dynamics of quantum walks and quantum 
lattice gas models Q has also been found. 

One of the key properties of quantum walks is their 
ability to evolve disentangled states into entangled ones 
and to efficiently generate entanglement in experimen- 
tally feasible systems Q. Controlled entanglement gen- 
eration has currently a place at the forefront of research, 
as it is a fundamental resource in quantum computation 
and cryptography and therefore a prc-requisite for the 
construction of reliable devices for quantum information 
processing Q . The role of entanglement in different pro- 
cesses in nature is also currently under investigation [lCf . 

Recently, an experimental realization of a linear quan- 
tum walk of a single neutral atom in a spin-dependent 
one-dimensional optical lattice has been reported, us- 
ing site-resolved fluorescence imaging to demonstrate the 
spatial coherence of the final wave function [ll|. Its 
conceptually straightforward generalization to higher- 
dimensions can however be experimentally demanding 
and would require significant technological efforts. Here, 
we propose a novel scheme for a two-dimensional quan- 
tum walk that significantly reduces the resources neces- 
sary for its realization. This is, in particular, true for 
an experimental setup similar to Ref. [Ill ]; our proposal 



able to mimic perfectly a specific spatial probability dis- 
tribution associated with a well-known two-dimensional 
quantum walk. Moreover, it is efficient in the genera- 
tion of spatial entanglement between the two orthogonal 
directions of the lattice on which the walker is moving. 

To introduce our notation, let us start by describing 
a discrete-time quantum walk (we will not deal with 
continous-timc quantum walks in this work) for the sim- 
pler one-dimensional model. In order to do that, we de- 
fine a two-dimensional Hilbert space, He (coin space), 
spanned by {|0), |1)}, and an infinite dimensional Hilbert 
space, Hw (walker space), spanned by {\i}}, with i as- 
suming all possible integer values. The state of the sys- 
tem is described as a vector in Hc®Hyj and the evolution 
of the system is given by a sequence of conditional shift 
and coin operations. We describe the conditional shift 
operation as 

S = |0> o (0|®5^ \i-l) w (i\+\l) c (l\®Yl HVwiil (1) 

iSZ iGZ 

If we consider the walker-component \i)w & s describing 
the quantized position of the walker on a line, with in- 
creasing numbers from left to right, the effect of S is to 
move the walker one step to the left (right) when the 
coin-component is in the state |0)c (|l)c)- In the orig- 
inal quantum walk Q the coin operation was chosen to 
be the Hadamard gate 
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can however be easily adapted to other schemes which are 
able to realize one-dimensional quantum walks. In the re- 
mainder of the paper, we demonstrate that our scheme is 



but it can be any unitary operation acting on He only. A 
single time step consists here of a coin operation followed 
by a shift operation. 

In the same way that classical random walks can be 
generalized to multi-dimensions, quantum walks do not 
have to be restricted to one dimension either and dif- 
ferent two-dimensional situations have been analyzed in 
literature [12|. However, the higher-dimensional walker 
space has required a similar increase in the dimension of 
the coin space in these schemes. In particular, for a two- 
dimensional walker, the coin system has been taken as 
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FIG. 1: (a) Spatial probability distribution after t — 50 steps 
of a two-dimensional Grover walk with the initial state of the 
coin as given in Eq. Q. Only the sites with even x and y 
are shown, as the probability is zero for all odd sites, (b) 
Sketch of the proposed scheme, in line with the experiment 
in Ref. [Tl|. A two-dimensional spin-dependent optical lattice 
can alternately be shifted into the two orthogonal directions. 



four-dimensional, which can equivalcntly be interpreted 
as two different two-level coins or a single four-level coin. 
For a simpler illustration, we consider the second case. 
The states of the computational basis of the coin \0)c, 
|1)C) |2)c an d |3)c then correspond to movements in 
the left-down, left-up, right-down and right-up direc- 
tions, respectively. The specific choice of the coin oper- 
ation distinguishes a number of different classes of two- 
dimensional quantum walks [l2j . 

Let us focus our attention on the two-dimensional 
quantum walk known as the Grover walk. This partic- 
ular walk has raised the interests of the scientific com- 
munity, as it can be used in order to implement the 
two-dimensional Grover search algorithm [13J. The cor- 
responding coin operation is given by 
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and it has been shown that for this specific coin operation 
the walker is always localized (i.e., the probability to find 
it at the origin is asymptotically larger than for t — > 00), 
except if the coin is in the particular initial state [Tij 
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The spatial probability distribution for this case can be 
obtained by tracing out the state of the coin and we show 
this distribution after t = 50 time steps in Fig. QJa). 

Unfortunately, these two-dimensional quantum walks 
are experimentally demanding. The use of a four-level 
quantum system or two distinct qubits (which requires 
however an entangling gate at each time step) as the 
coin complicates the realization and only a few feasi- 
ble physical implementations have been proposed so far 
(see, for instance, Ref. |f|). In order to reduce the ex- 
perimental challenges, we suggest considering a different 
two-dimensional quantum walk, in which the coin is a sin- 
gle qubit, and the movement on the x and y directions 



are alternate [according to the sequence: coin operation 
- movement oni- coin operation - movement on y, see 
Fig. Hfb)]. For this, we start with the initial state of the 
coin 
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(5) 



as in the original symmetric one-dimensional quantum 
walk, and consider the Hadamard gate as the coin oper- 
ation. A time step consists here of a sequence of the two 
Hadamard operations and the two movement on the x 
and y directions. Surprisingly, if we analyze the spatial 
probability distribution of this two-dimensional quantum 
walk, we find exactly the same result as for the non- 
localized Grover walk described above [see Fig.QJa)]. 

To prove formally this equivalence, in what follows we 
will show how the coefficients of the Grover walk in the 
non-localized case can be mapped to the coefficients of 
the alternate walk state with the initial condition of the 
coin given by Eq. (|5|). For this, we define the basis states 
of the Grover walk as {\x,y,c}}, where x, y G Z denote 
the position along the x and y directions, respectively, 
and c G {0, 1, 2, 3} describes the state of the coin. At the 
same time, we decompose the state of our alternate walk 
over the basis {|a;, y, c')}, where d G {0, 1} is the state of 
the single-qubit coin. The coefficients in the decomposi- 
tion of the Grover walk and of the alternate walk are then 
defined as a X:V:C (i) and /3 X)I/)C '(i), respectively. It is easy 
to note that, for the initial states under consideration, 
the a x ,2/,c(D's are real numbers, while the f3 x _ y _ c i(tys are 
complex numbers. We will now show that 

P*,v,o(t) = (-D*e l7r/4 K, y ,oW + ia x , y ,2(t)], (6) 
Px,vM = (-l)*e i " / *[-a a ,, v ,i(t) + *a*,v,3(i)]. (7) 

For this, let us first demonstrate that the amplitudes sat- 
isfy the properties 

a x -i. y .o(t) + a x -i. y .i(t) + a x+ i iVi 2(t) + a x+ i. y . 3 (t) = 0, 

(8) 

Cix,y-l,o{t) + a x , y -l t 2(t) + a&, y +i,i(t) + a x ,y+i,a{t) = 0, 

(9) 

if the walk starts at the origin with the initial state cor- 
responding to 



«o.o,o(0) = 1/2, ao,o,i(0) = -1/2, 
00,0,2(0) = -1/2, 00,0,3(0) - 1/2 



(10) 



and all other a's at t = being zero. 

Our proof works by induction on t. It is easy to verify, 
by means of a direct calculation, that Eqs. © and © 
are satisfied at t = 0. For this, we only need to evaluate 
them for the values (x,y) equal to (1,0), (—1,0), (0,1) 
and (0,-1), as for all other values the a's arc initially 
zero. Then, we assume that both Eqs. (jHJ) and © are 
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true for any couple of values (x, y) at some time t, and 
we need to prove that they hold at time t + 1 . Starting 
with the left hand side of Eq. © , we have 
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= ^2 a Xl y-l t 2i(t) + a x,y+l,2i+l(t), (11) 

which is identically zero because we have assumed Eq. © 
true at time t. Here, G^ = 0, 1, 2, 3) is the element 
of the matrix G corresponding to |«)c(j|- It is straight- 
forward to proceed in the same way to prove Eq. §§§ at 
time t + 1, assuming that Eq. ((SJ) is true at time t. 

We can now prove the relations in Eqs. (O and (J7J for 
the initial conditions for the alternate walk given by 
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Aai(0)='/A (12) 



and all other /3's being zero at t = and Eq. (|TU1) for 
the Grover walk. Again, we proceed by induction in t. 
At t = 0, all the amplitudes a's and /3's are zero outside 
of the origin (a;, y) = (0, 0) and using Eqs. (flQ|) and (fT2|) 
one can straightforwardly verify that Eqs. (j6|) and (0 arc 
fulfilled at the origin. We then assume that Eqs. (JB]) and 
([7]) are true [for any couple of values (x, y)] at some time 
t, and we prove that they hold for the time t+1. For our 
alternate walk this leads to 

Px,y,i{t + 1) = ^{Px+l,y+(-lY + Px+l,y+{-Xy,\{t) 

+ (-l) < [j9x-i,v + (- 1 )i,o(t)-i9x-i,„+(-i)M(*)]} ) (13) 
with i = 0, 1, while, for the Grover walk, we find 

3 

a x , v ,k(t + 1) = E G fei a x+(-i) m ,a+(-i)" ( 14 ) 

3=0 

with fc = 0,1,2,3, to = Int[fc/2] and n = mod[2,A;]. 
Starting from Eq. (|13l) with i = and assuming Eq. © 
true at time t, we can therefore write 

Px,y,o(t + l) 

1 1 1 

= 2(- 1 ) vi {E(- 1 ) <a -+ i '»+ i .«(*)+E a «- i .»+ i .«(*) 
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(15) 



Now, by using Eq. 
to 

i8x,»,o(t + l) 



straightforward calculations lead 
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= (-l)* +1 e ? T ^ [GoiOx+i.v+x.jW + iGaj-a^-i^+i,^*) 

= (-l) t+1 e^/ 4 [a a , tf , (« + 1) + ia a , tf , a (t + 1)] , (16) 

which completes the proof for Eq. (JBJ. An analogous 
analysis allows us to prove the partner Eq. ([7]), showing 
that the alternate walk and the Grover walk generate the 
same spatial density distribution for the aformentioncd 
initial conditions. 

As already stated, one interesting aspect of quantum 
walks is their ability to generate entanglement. Even if 
the spatial probability distributions are the same in the 
two different walks investigated in this paper, the final 
state after a fixed number of steps is clearly different. 
It is therefore interesting to make a comparison between 
them in terms of entanglement generation. A well inves- 
tigated feature in one- as well as two-dimensional quan- 
tum walks is the generation of entanglement between the 
state of the coin and t he p osition of the walker (coin- 
position entanglement) [15|. Since the evolution of the 
whole system under the action of the walk is unitary, the 
total state remains pure, and it is possible to use the 
von Neumann entropy of the reduced density matrices 
to measure the coin-position entanglement. For this, we 
need to evaluate S(pc) = S(pw), where pc (pw) is the 
reduced density matrix of the coin (walker). To make 
a comparison between the two different dynamical walks 
described above, we first calculate S(pc) for the Grover 
walk, with the initial conditions as in Eq. (flQ|) . The re- 
sults for t — - 20 steps are shown in Fig. [2ja). One can 
clearly see that the entropy is quite high throughout the 
process and close to its maximum, which in this case is 
S = log 2 4 = 2, because the Hilbert space of the coin is 
four-dimensional. The coin-position entanglement in the 
alternate walk, with the initial conditions as in Eq. (|12l) , 
is also shown in the same figure. Similarly to the Grover 
walk, the entropy quickly reaches values close to its maxi- 
mum (the fluctuations are less significant here). However, 
due to the fact that the Hilbert space of the coin is only 
two-dimensional, this corresponds to S = log 2 2 = 1. 

Let us now consider a different kind of entanglement. 
Suppose that we want to entangle the orthogonal direc- 
tions of the lattice on which the walker is moving (we will 
denote this as x-y spatial entanglement). In order to do 
so, we need to trace out the degree of freedom embodied 
by the coin. This can be an advantage of our alternate 
walk in terms of the amount of entanglement generated, 
due to the smaller dimension of the coin space. How- 
ever, since the state after tracing out the coin is mixed, 
the von Neumann entropy can no longer be used as a 
measure of entanglement and we will, for this reason, use 
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FIG. 2: (a) Coin-position entanglement against the number 
of steps t for the Grover walk (blue circles) with the initial 
conditions as in Eq. (|10p and for our alternate walk (purple 
squares) with the initial conditions as in Eq. (|12[) . (b) En- 
tanglement between x and y position of the particle for the 
Grover walk and our alternate walk, with the same conditions 
as in panel (a). 



we have sketched a possible physical implementation of 
the scheme, which is in lines with the current experimen- 
tal technology. A deeper investigation of two-dimensional 
quantum walks with singlc-qubit coins will be an exciting 
and interesting extension of the work presented here, due 
to the importance of quantum walks in the implementa- 
tion of quantum algorithms (for instance, an extension to 
a more general coin operation, along the lines of Ref. , 
is currently under investigation.). We believe that our re- 
sults contribute to this task in a significant way. 
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sions. This work was supported by Science Foundation 
Ireland under grant no. 05/IN/I852 and 05/IN/I852 NS. 
C.D.F. is supported by the Irish Research Council for 
Science, Engineering and Technology. 



the negativity N of the partial transpose, in its gener- 
alization for hig her-dimensional systems (so as to have 
< N < 1) [l6[. The results for both walks for a num- 
ber of steps up to t = 20 are shown in Fig.[2][b), and one 
can clearly see that for the alternate walk the x-y spatial 
entanglement is, in fact, always larger. 

Let us finally briefly discuss a possible experimental 
implementation of our proposed walk, following Ref. . 
In the experiment of Karski et ah, a single cesium atom is 
able to walk in a one-dimensional optical lattice in posi- 
tion space. The conditional shift operator is realized by 
continuously controlling the trap polarization in a way 
so as to move the spin state |0) (|1)) to the right (left). 
Coin operations are obtained by proper laser pulses and 
the atom distribution is probed by fluorescence imaging. 
The setting can be adapted to our walk by keeping the 
same coin operation and allowing the optical lattice (in 
this two-dimensional one) to alternately shift into 

the two orthogonal directions, as sketched in Fig. QJb). 
In terms of experimental challenges, this scheme is there- 
fore more feasible and could in principle be realized with 
the state-of-the-art current technology. In fact, it just ex- 
ploits a straightforward extension of the existing setup 
and is thus less demanding than performing a coin oper- 
ation on four different internal states of the particle and 
moving them into four possible directions simultaneously 
(as required for an implementation of the original Grover 
walk) . 

In summary, we have shown that a two-dimensional 
quantum walk with a four-dimensional coin space can 
be equivalent to two one-dimensional walks in alternat- 
ing directions at alternating time steps, for which only 
a two-dimensional coin space is required. In particular, 
we have formally shown that such a walk is able to per- 
fectly mimic the spatial probability distribution of the 
well-known Grover walk. We have analyzed the behav- 
ior of two different kinds of entanglement for both walks 
and found that the presented alternate walk is more ef- 
ficient at generating spatial x-y entanglement. Finally, 



[1] M. Barber and B. W. Ninham, Random and Restricted 
Walks: Theory and Applications (Gordon and Breach, 
New York, 1970); B. Malkiel, A Random Walk Down 
Wall Street (W. W. Norton, New York, 1973); H. C. 
Berg, Random Walks in Biology (Princeton University 
Press, Princeton, NJ, 1993). 

[2] Y. Aharonov, L. Davidovich, and N. Zagury, Phys. Rev. 
A 48, 1687 (1993). 

[3] J. Kempe, Contemp. Phys. 44, 307 (2003). 

[4] A. M. Childs, Phys. Rev. Lett. 102, 180501 (2009). 

[5] B. C. Travaglione and G. J. Milburn, Phys. Rev. A 65, 
032310 (2002); W. Diir et al, Phys. Rev. A 66 052319 
(2002); Z. Zhao et al, |arXiv:quant-ph/0212149| (2002); 

B. C. Sanders et ah, Phys. Rev. A 67, 042305 (2003); H. 
Jeong, M. Paternostro, and M. S. Kim, Phys. Rev. A 69, 
012310 (2004). 

[6] K. Eckert et al, Phys. Rev. A 72, 012327 (2005). 

[7] R. P. Feynman and A. Hibbs, Quantum Mechanics and 

Path Integrals (McGraw-Hill, New York, 1965). 
[8] S. E. Venegas-Andraca and S. Bose, larXiv:0901.3946l 

(2009); S. K. Goya and C. M. Chandrashekar, J. Phys. 

A 43, 235303 (2010). 
[9] R. Horodecki et al., Rev. Mod. Phys. 81, 865 (2009). 
[10] G. S. Engel et al, Nature 446, 782 (2007); H. Lee, Y.- 

C. Cheng, G. R. Fleming, Science 316, 1462 (2007); M. 
Mohseni et al, J. Chem. Phys. 129, 174106 (2008); M. B. 
Plenio and S. F. Huelga, New J. Phys. 10, 113019 (2008); 
J. Cai, G. G. Guerreschi, and H.-J. Briegel, Phys. Rev. 
Lett. 104, 220502 (2010). 

[11] M. Karski et al, Science 325, 174 (2009). 

[12] T. D. Mackay et al, J. Phys. A 35, 2745 (2002); A. C. 

Oliveira, R. Portugal, and R. Donangelo, Phys. Rev. A 

74, 012312 (2006). 
[13] N. Shenvi, J. Kempe, and K. B. Whaley, Phys. Rev. 

A 67, 052307 (2003); A. Ambainis, J. Kempe, and A. 

Rivosh, in Proc. 16th ACM-SIAM SODA, Vancouver 

(SIAM, Philadelphia, USA, 2005), p. 1099; A. Tulsi, 

Phys. Rev. A 78, 012310 (2008). 
[14] N. Inui, Y. Konishi, and N. Konno, Phys. Rev. A 69, 

052323 (2004); M. Stefanak, T. Kiss, and I. Jex, Phys. 

Rev. A 78, 032306 (2008). 
[15] I. Carneiro et al, New J. Phys. 7, 156 (2005); G. Abal 



5 



et al, Phys. Rev. A 73, 042302 (2006); M. Annabestani, [16] S. Lee et al., Phys. Rev. A 68, 062304 (2003). 

M. R. Abolhasani, and G. Abal, J. Phys. A 43, 075301 [17] K. Watabe et al., Phys. Rev. A 77, 062331 (2008). 

(2010). 



